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Origin Of The Cosmic Ray Positrons Observed Near Earth- 
Meson Decay Or Dark Matter Decay? 
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ABSTRACT 

We show that the flux of the high energy cosmic ray positrons observed near 
Earth is that expected from the decay of mesons produced by the primary cosmic 
rays (CRs) in the local interstellar medium. 

Subject headings: cosmic rays, dark matter 


1. Introduction 

A wide variety of evidence points to the existence of dark matter in the universe which 
cannot be seen directly, but which can be detected by its gravity. Precise measurements of 
the angular power spectrum of the cosmic microwave background radiation analyzed with 
the standard model of cosmology indicate that nearly 27% of the total mass-energy density 
of the universe is residing in this dark matter (Bennet et al. 2013; Ade et al. 2014) whose 
existence was first discovered in clusters of galaxies (Zwicky 1933,1937) and later in galaxies 
(Rubin 1970). The nature and origin of this dark matter (DM) are still unknown. 

If the dark matter is made of weakly interacting massive particles (WIMPs) relics from 
the big bang, their annihilation cross section satisfies < cx w v >~ 2 x 10 -26 cm 3 /s (Jungman 
et al. 1996), which is typical of weakly interacting ~TeV-mass particles. Annihilation/decay 
of such particles can produce stable particles such as electrons, positrons, protons, antipro¬ 
tons, neutrinos and gamma rays, whose energy cannot exceed the mass of the DM particles. 
Hence, in recent years, bumps with a sharp cutoff due to the decay and/or annihilation 
of WIMPs were looked for in high precision measurements of the energy spectrum of such 
cosmic ray (CR) particles. A bump around E=620 GeV with a sharp cutoff around E=800 
GeV in the combined flux <f> e ± of e 1 * 1 cosmic rays (CRs) was reported by the Advanced Thin 
Ionization Calorimeter (ATIC) balloon experiment (Chang et al. 2008). It was interpreted 
as a possible dark matter signal (e.g., Bergstrom et al. 2008; Cholis et al. 2009; Circlli et 
al. 2009; Arkani-Hamed et al. 2009; Pohl 2009). However, such a peak was not confirmed 
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by measurements with the ground based Cherenkov Telescopes of High Energy Stereoscopic 
System (H.E.S.S.) (Aharonian et al. 2009), and later with more precise measurements with 
the Large Area Telescope (LAT) aboard the Fermi satellite (Ackermann et al. 2010) and 
with the Alpha Magnetic Spectrometer (AMS) on the International Space Station (Aguilar 
et al. 2014). But, at the same time, the energy spectrum of e 1 * 1 CRs measured with H.E.S.S. 
(Aharonian et al. 2008) showed a sharp break near 1 TeV whose origin is still not clear. 

Shortly after the ATIC report (Chang et al. 2008), an unexpected increase of the CR 
positron fraction <L(e + )/(<L(e + ) + <L(e - )) as a function of energy between ~10 GeV and ~100 
GeV, observed with the Payload for Antimatter Matter Exploration and Light-nuclei Astro¬ 
physics (PAMELA) satellite was reported (Adriani et al. 2009,2010,2013). Unlike the ATIC 
bump, the rise with energy of the positron fraction has been confirmed by measurements 
of the separate CR ray fluxes of e 1 * 1 with Fermi-LAT (Ackermann et al. 2012), and with 
AMS (Aguilar et al. 2013,2014; Accardo et al. 2014) with much larger statistics and smaller 
systematic errors. In the high precision AMS data, this rise continues beyond 100 GeV but 
with a decreasing rate that appears to level off around 275 ± 32 GeV (Accardo et al. 2014, 
Aguilar et al. 2014, Pohl 2014). 

The measured flux of CR positrons and the increasing positron fraction were in stark 
contrast with those predicted for secondary production of positrons in the interstellar medium 
(ISM) by primary CRs from detailed calculations with the elaborate GALPROP code (e.g., 
Moskalenko and Strong 1998a,b). This sparked the publication of several alternative expla¬ 
nations to the origin of the excess of high energy CR e + flux, such as dark matter (e.g., 
Bergstrom et al. 2008,2009,2013; Bergstrom 2009; for a recent review see, e.g., Ibarra et al. 
2013 and references therein), positron emission from a few nearby pulsars and/or supernova 
remnants (e.g., Hooper et al. 2009; Shaviv et al. 2009), secondary production of positrons 
in collisions of primary CRs in/near their sources (Dado & Dar 2010) such as the highly rel¬ 
ativistic jets launched mainly by supernovae Type Ic and acreting stellar mass and massive 
black holes (e.g., Dar and De Rujula 2004, 2008 and references therein) and/or the spherical 
ejecta in supernova remnants (SNRs) (Blasi 2009; Stawarz et al. 2010). Other authors (e.g., 
Blum et al. 2013; Cowsik et al. 2014) maintained that secondary production of positrons in 
the ISM by the primary CRs can explain the observed flux of CR positrons although their 
calculations involved incorrect assumptions concerning the energy-loss of positrons in the 
ISM. 


The reported isotropy of the flux of high energy CR positrons observed with AMS 
(Accardo et al. 2014), however, is in tension with the assumption that the observed flux of 
e + CRs was produced by a few relatively nearby pulsars (e.g., Linden and Profumo 2013) 
or supernova remnants. At the same time, the high precision data of AMS02 (Aguilar et 
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al. 2014) on the combined e^ flux up to 1 TeV and the separate e~ and e + fluxes, and the 
positron fraction (Accardo et al. 2014, Aguilar et al. 2014) that appear to flatten above 200 
GeV, do not show yet any convincing evidence for a cutoff/sharp decline, which begins below 
~ 700 GeV that could be associated with a decay/annihilation of dark matter, particles. 

In this paper, we re-examine the possibility that the main origin of the high energy CR 
positrons observed near Earth is secondary production in hadronic collisions of the primary 
cosmic ray protons and nuclei in the ISM (Blum et al. 2013; Cowsik et al. 2014). We 
focus on the high energy (E> 10 GeV) behavior of the flux of CR e +, s where it may be 
dominated by the decay/annihilation of dark matter particles. We consider in detail, the 
energy-loss of e + ’s above 10 GeV where it is dominated by synchrotron radiation, inverse 
Compton scattering, triple pair production, and escape from the Galaxy by diffusion in 
its turbulent magnetic fields. (Other energy loss mechanisms that are important only at 
energy well below 10 GeV, such as Coulomb scattering, ionization and bremsstrahlung, as 
well as threshold effects, geomagnetic shielding and solar modulation, are not considered in 
detail, but are included for completeness through a best fit phenomenological depletion factor 
D e ±(E) = 1 — exp(—(E/V) Q ), which depends on the time of the measurements and on the 
clecctric charge (±) but does not affect the behavior at E>10 GeV.) Using only standard 
particle physics, the steady state approximation for CR propagation, and the measured 
properties of the primary CRs and the local Galactic environment, we show that the spectrum 
of high energy CR e + ’s measured by AMS (Aguilar et al. 2014) with high statistics and small 
systematic errors is consistent with that expected from secondary production of positrons in 
hadronic interactions of the primary CRs in the local ISM. In particular, the hardening of 
their spectrum with increasing energy between ~ 30 and 150 GeV and its softening beyond 
250 GeV, as observed with AMS02 (Aguilar et al. 2014), are expected from the transition of 
their inverse Compton scattering of diffuse Galactic light from the Thomson regime to the 
Klcin-Nishina regime (Schlickeiser & Ruppel 2010). 

Moreover, the break near 1 TeV in the combined e 1 * 1 flux, suggested by the data from 
H.E.S.S. (Aharonian et al. 2008) may be the maximum re-acceleration energy of electrons in 
the ISM due to their radiative energy-loss and escape from the Galaxy by diffusion. Indeed a 
break was reproduced in numerical calculations (e.g., Stawarz et al. 2010), which have used 
the elaborate GALPROP numerical code (Moskalenko & Strong 1998a,b). Other standard 
astroparticle physics interpretations rather than a dark matter annihilation/decay signal are 
shortly discussed. 
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2. CR Production of positrons in the local ISM: 

The flux of primary Cosmic ray nucleons (free protons and nucleons bound in atomic 
nuclei) in the energy range between several GeV and PeV per nucleon observed near Earth 
is well described by (e.g., Olive et al. 2014), 

4>p(E) « 1.8 (E/GeV) -/3 fu, (1) 

where /3 = 2.7 and fu = l/(GeVcm 2 s sr) is the flux unit. The contribution of nucleons 
bound in atomic nuclei (A,Z) to the cosmic ray flux at the same energy per nucleon is less 
than 5% and comes mainly from light nuclei whose inelastic cross section per nucleon is 
~ <jp a/A ~ cr p p, i.e., roughly that of free protons, while their effective rigidity (R=AE/Z) is 
twice than that of free protons. Hence, due to Feynman scaling (1969), secondary production 
of high energy positrons through the decay of n and K mesons, which are produced in 
hadronic collisions of cosmic ray nucleons with the baryons in the local ISM where the mean 
baryon density is n ism « 0.9 cm -3 (Kalberla & Dedes, 2008) inject into the CR halo a flux of 
e + ’s per unit volume roughly at a rate 

J e + ~ Ke+ ^in(pp) C llj sm 4> p (E) . (2) 

where K e + is a constant that depends only on the power-law index of the product cr in (pp) <3> p (E) 
bnt not on E itself. In a steady state, the flux <3> e +(E) of secondary e +, s produced in the 
ISM satisfies 

A[ b (E)4v(E)] = J« + (E) (3) 

where b(E) = dE/dt is the loss rate of e + energy by radiation (rad) and by escape (esc) from 
the Galaxy by diffusion through its turbulent magnetic field. The solution of Eq. (3) for 
$p oc E -/3 is 

4> e + (E) « K e + cr in (pp) n ism c r e $ p /(/3j - 1) • (4) 

where a in (pp) ~ 30 (E/GeV) 0 ' 06 mb, K e + « 7 x 10 -3 for /3j « f3 — 0.06 = 2.64, and r e = E/(dE/dt) 
is the mean life-time of positrons in the ISM due to their escape (esc) from the Galaxy by dif¬ 
fusion and radiative (rad) energy losses (inverse Compton scattering of background photons 
and synchrotron radiation). Hence, the expected flux of high energy CR e +, s near Earth is 
simply 

4> e+ (E) « 6.22 x 10 -18 E -2 ' 64 r e (E) fu. (5) 

where the combined radiative losses and escape by diffusion yields 


l/^e 1/^esc T 1/'Tr ac i . 


(6) 



5 


At high energy, the rigidities of e +, s and protons become practically equal. Consequently, 
their escape times by diffusion become practically equal. For random Galactic magnetic 
fields with a Kolmogorov spectrum 

r esc « (7.5 ± 1.5) x 10 14 (E/GeV)' 1/3 s (7) 

where the normalization has been adjusted to the value obtained by Lipari (2014) from a 
leaky box model analysis of the flux ratio 10 Be/ 9 Be measured with the Cosmic Ray Isotope 
Spectrometer (GRIS) in the energy range 70-145 MeV/nucleon (Yanasak et al. 2001). At 
the CR ankle near E w 4 x 10 9 GeV, this normalization yields r esc = (15 ± 3) ky, consistent 
with an expected free escape time H/c=5 ± 1 kpc/c=(16.3 ± 3.2) ky, where H is the typical 
scale height of the Galactic cosmic ray halo perpendicular to the Galactic disk (e.g., Trotta 
et al. 2011). 


Synchrotron emission and inverse Compton scattering of background photons in the 
Thomson regime yield 

T rad « ~ 0.67 x 10 16 (E/GeV)- 1 s (8) 

4 cr T c U E 

where a T ~0.66 barn is the Thomson cross section and U ~ 1.47eV/cm 3 is the energy den¬ 
sity of the background radiations plus the magnetic field in the local ISM. Roughly, the 
magnetic field B ~ 4/iG in the local ISM contributes B 2 /8 7r ~ 0.40eV/cm 3 to U, the dif¬ 
fuse Galactic light (DGL) contributes ~ 0.39eV/cm 3 , the far infrared (FIR) light contributes 
~ 0.40eV/cm 3 , and the cosmic microwave background (CMB) contributes 0.26eV/cm 3 (e.g., 
Porter et al. 2006; Schlickeiser & Ruppel 2010). 


In the Klein-Nishina regime, (w = 4E e e 7 /(m c c 2 ) 2 >l), the energy loss-rate by ICS is 
given by (e.g., Blumenthal & Gould 1970) 

dE 3(t t cU 7 (m e c 2 ) 2 
dt 8 e 2 


■[lnw — 11/6] 


( 9 ) 


Consequently, as can be seen from Eq. (3), in a steady state, an injected power-law spectrum 
of electrons becomes softer in the Thomson regime, E _/3j —> E _/3j_1 , whereas in the Klein- 
Nishina regime it becomes harder, E~ /5 j — > E” /3j+1 . 


In the case of a gray body radiation of temperature T, the transition from the Thomson 
regime to the Klein-Nishina regime can be interpolated through replacing U(E) by 


Ukn(E) 

where (Schlickeiser & Ruppel 2010) 


U(E) E| n /(E 2 + E 2 N ) 


2\2 


E 


ICN 


0.27 (m e c 2 ) 

k¥ 


( 10 ) 


( 11 ) 
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with k being the Boltzman constant. For the local DGL, Tfs 5700°K and E K n ~ 140 GeV. 
For the Galactic FIR radiation and the CMB, Ekn is well above 10 TeV. 

Since the estimated values of the local energy densities of the DGL and FIR, as well 
as the local B suffer from uncertainties, we have best fitted their values within the er¬ 
rors of their estimated values. At 100 GeV, they yield for instance, r esc = 1.62 x 10 14 s, 
T r ad ~ 7.2 x 10 13 s, and consequently r e ~ 0.50 x 10 14 s ~ 1.6 MY. 

The flux of high energy e + CRs measured with AMS02 (Aguilar et al. 2014) and the 
flux expected from CR production of e +, s in the ISM as given by Eq. (5) multiplied by 
the low-energy best fit phenomenological depletion factor D e +(E) = 1 — exp(—(E/V) Q ) with 
V=3.90 GeV and a = 1.78 are compared in Fig. 1. 


3. Origin of the CR e flux 

Most of the Galactic accelerators of high energy cosmic rays (supernova remnants, pul¬ 
sars, gamma ray bursts, and the massive black hole near the Galactic center) are located in 
the Galactic disk within ~ 4 kpc from its center (e.g., Case & Bhattacharya 1998; Yusifov, 
& I. Kucuk 2004; Lorimcr et al. 2006). Presumably, the Fermi acceleration of high energy 
electrons, protons and nuclei is simultaneous and imparts to them the same Lorentz factor 
distribution (e.g. Dar & De Rujula 2008), i.e., a source distribution J e - oc J p . These primary 
CRs can reach Earth from their common sources by diffusion (dif) through the turbulent 
Galactic magnetic fields. High energy CR protons and electrons have practically the same 
rigidity (except for the sign). Consequently, without energy losses, they would have reached 
Earth from their common sources after the same mean diffusion time, and their fluxes in the 
local ISM then would have satisfied 0 e -(E) oc <L P (E). This relation can hold as long as their 
energy loss by radiation is quite small during their diffusion time from source to Earth. 

The radiative energy loss of e _ CRs in the ISM satisfies dE/dt = —bE 2 whose solution 
is E = E 0 [1 — bEtdif] where E is the observed energy of the e _ CR near Earth and E 0 is its 
injection energy. It implies that most of the sources that are located beyond ~ 4 kpc cannot 
contribute much to the observed local flux of CR e - ’, unless their CRs are re-accelerated in 
the ISM or transported ballistically over long Galactic distances by highly relativistic jets, 
such as those emitted in supernova explosions of Type Ic, most of which do not point in our 
direction (e.g., Dar & Plaga 1999, Dar & De Rujula 2008). 

If the escape time of CR protons and electrons from their sources into the ISM is much 
shorter than their energy loss rate there, then the Fermi accelerated high energy electrons 
are injected into the ISM at a rate J e oc < h p /r esc (E). Hence, their steady state flux near Earth 
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satisfies 

$ e - = D e - (E) A e - E~ p T e /r esc + 0.74 4> e + , (12) 

The first term on the right hand side (RHS) of Eq. (12) is due to primary CR e _ ’s, presumably 
re-accelerated in the ISM or transported ballistically by highly relativistic jets from all distant 
sources together with the CR protons (A e - is an unknown flux normalization constant, 
which is treated as a free parameter.) The second term on the RHS is the contribution 
from hadronic production of high energy e _ ’s in pp collisions, in the local ISM assuming 
an inclusive production ratio identical to that of muons, e _ /e + ~ yU~//i + = 0.74, which was 
measured for atmospheric muons and in accelerator experiments (Haino et al. 2004; Archard 
et al. 2004; Adamson et al. 2007; Khachatryan et al. 2010; Agafonova et al. 2010). Eq. (12) 
can also be written as 


E 3 $ e - 


D e - (E) A e - E 3 ~ 0 

(1 T ^esc/^rad) 


+ 0.74E 3 $ e+ , 


(13) 


Assuming that the mean radiation field in the CR halo can be represented by that in the 
solar neighborhood, the predicted flux of CR e~’s as given by Eq. (13) for the best fit 
values A e - = 0.27GeV 3 fu, V=2.25 GeV and a = 1.46 is compared in Fig. 2 to the CR 
e~ flux measured with AMS02 (Aguilar et al. 2014). The agreement is quite satisfactory 
(x 2 /df = 0.85). 


4. The e + fraction and the combined e 11 " flux 

The positron fraction <3> e +/<I>e as a function of energy obtained from our predicted e 1 * 1 
fluxes, which are plotted in Figs. 1 and 2, and the positron fraction measured with AMS02 
(Accardo et al. 2014) near Earth are compared in Fig. 3 (y 2 /df = 0.35). The combined e 1 * 1 flux 
4> e = 4> e + + <3> e - measured near Earth with AMS02 (Aguilar et al. 2014) and the expected 
flux obtained from standard astroparticle physics are compared in Fig. 4 (y 2 /df = 1.27). As 
can be seen from Figs. 1-4, the agreement between the near Earth e 1 * 1 CR fluxes measured 
with AMS02 and those expected from standard astroparticle physics is quite good. 


5. Dark matter annihilation signal in AMS02 ? 

The AMS collaboration presented an impressive ’minimal model’ fit (y 2 /df = 0.63) to 
their updated data (Accardo et al. 2014) on the CR positron fraction, which is shown in 
Fig. 5. In this model the e 1 * 1 CR fluxes were parametrized as the sum of a power law 
spectrum and a common source term with an exponential cutoff, which may represent a 



dark matter decay/annihilation contribution, 


$ e+ = Ce+ E _7e+ + C s e _E ^ Es , 
<h e - =Ce- E“ 7e - +C s e- E / Es , 


(14) 

(15) 


with E in GeV. However, the positron fraction < h e +/( < h e - + ^e+) is invariant under division 
of both fluxes by an arbitrary E-dependent function. In particular, by dividing the minimal 
model fluxes by C e + E _7 e+, a very good best fit (y 2 /d/ = 36.4/58) to the positron fraction 
in the energy range 1 to 500 GeV yielded (Accardo et al. 2014) C e +/C e - = 0.091 ±0.001, 
7e+ — 7e- = 0.56 ± 0.03, C s /C e - = 0.0061 ± 0.0009, 7 S — %- = —0.72 ± 0.04 and a cutoff pa¬ 
rameter E s = 0.54 ± 0.17 TeV. These best fit parameters specify the e 1 * 1 fluxes of the minimal 
model only up to an unknown common power-law factor C e +E _7 e+. 


A major feature of the minimal model is the sharp decline of the positron fraction above 
~500 GeV. In order to test the AMS minimal model, we have best fitted $6+ and <f> e - as 
given by Eqs. (13) and (14), respectively, to the separate e + and e _ fluxes measured with 
AMS02 (Aguilar et al. 2014). The best fit parameters were constrained to reproduce the 
above relations between the parameters of the minimal model. Unlike the minimal model 
impressive fit to the positron fraction published by the AMS collaboration, (Accardo et 
al. 2014) as shown in Fig. 5, the best fits to the separate e + and e _ fluxes, and to their sum, 
are rather poor, as shown in Figs. 6-8. That is true in particular for E< 10 GeV and E> 300 
GeV. 


One may argue, however, that the minimal model represents well the source spectra of 
e ± CRs, which are later modified by propagation effects in/near source, in the ISM and in 
the heliosphere. But, because solar modulation, escape by diffusion in the turbulent Galactic 
magnetic fields, and energy-loss through ionization, bremsstrhalung, synchrotron radiation, 
inverse Compton scattering and pair production, which strongly affect the observed spectra 
of CRs near Earth, are independent of the sign of charge, they modify the injection spectra 
of e ± by the same energy-loss factor E/(dE/dt) = b(E). However, the spatial distribution 
of Galactic dark matter is quite different from that of the more conventional sources of 
primary e - and secondary e^ CRs. This makes rather unlikely the possibility that the 
strong modifications of the source spectra of both high energy e + and e _ CRs observed near 
Earth are nearly identical. 


6. CR e 1 * 1 knee near TeV? 

The energy spectrum of high energy e 1 * 1 CRs measured with H.E.S.S. (Aharonian et 
al. 2008,2009) suggests a sharp break in the combined e 1 * 1 spectrum near E~TeV. This is 
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shown in Fig. 10 where we plotted a smooth cutoff power-law fit to the combined e 1 * 1 flux 
measured with AMS02 (Aguilar et al. 2014) and with H.E.S.S (Aharonian et al 2008, 2009) 
after normalizing the H.E.S.S data within their reported systematic errors to match the more 
precise AMS02 sub-TeV data. This fit yield a break energy around 1 TeV. 

Fermi acceleration of electrons is cut-off by synchrotron emission when the energy loss 
rate by synchrotron radiation exceeds the energy gain by magnetic deflections. This syn¬ 
chrotron cutoff energy is given roughly by E « a/ 6 e/<r T B m c c 2 . However, the typical mag¬ 
netic field B in supernova remnants (SNRs), which are widely accepted as the main source 
of high energy Galactic CR electrons, is well below 100 pG. The corresponding synchrotron 
cutoff in the energy spectrum of CR e’s accelerated in SNRs is expected at > 1 PeV, well 
above that suggested by the H.E.S.S. data. 

The observed flux of primary e _ CRs is also cut-off when their radiative life becomes 
shorter than their travel time by diffusion from source to Earth as discussed in Section 3. 
However, for the main Galactic sources, this cut-off is much below the H.E.S.S. break/cut-off. 

The H.E.S.S. cutoff could be a dark matter signal. However, the data of Fermi-LAT 
(Ackermann et al. 2010) and AMS02 (Aguilar et al. 2014) published so far, do not show any 
hint for such a break below TeV. Moreover, such a break, if real, still could be explained by 
standard astroparticlc physics rather than as a dark matter signal. The origin of a break 
could be, e.g., 

(A) A re-acceleration cutoff when the re-acceleration time in the Galactic ISM exceeds the 
electrons’ life time due to their radiative energy losses and escape from the Galaxy by 
diffusion. 

(B) CR e _ ’’knee” at E knee (e _ ) = (m e /m p )E knee (p) ~ 1 TeV in the spectrum of e _ CRs, which 
are Fermi accelerated together with protons and nuclei by the highly relativistic jets launched 
in SNe Ic to the same Lorentz factor distribution (e.g., Dar and De Rujula 2008). (E knee (p) ~ 2 
PeV is the CR ’’knee” in the flux per nucleon of CR nuclei). 

(C) A cutoff when the energy deposition rate of high energy electrons in the ISM by 
highly relativistic jets cannot compensate anymore their energy loss by radiation and escape 
from the Galaxy by diffusion. 

For other possible standard astroparticle physics origins of the H.E.S.S break see, e.g., 
Stawarz et al. (2010) and Cowsik et al. (2014). 
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7. Conclusions 

Figs. 1-4 demonstrate that standard astroparticle physics, such as positron production 
in the ISM in hadronic collisions of the primary CRs with gas in the local ISM and electron 
acceleration in a variety of Galactic sources, can explain the observed fluxes of CR electrons 
and positrons near Earth, which were measured with AMS02 (Accardo et ah 2014; Aguilar et 
al. 2014) without invoking a contribution from annihilation/decay of dark matter particles. 
In particular, the hardening of the positron flux between 20 and 200 GeV can be due to the 
Klcin-Nishina supression of the energy loss of positrons by inverse Compton scattering from 
the diffuse Galactic light in the local ISM. 

Neither the current published data on the spectra of the separate and combined e 1 * 1 CRs 
near Earth, nor the minimal model fit by the AMS collaboration to their measured positron 
fraction (Accardo et ah 2014) provide compelling evidence or hints for a contribution from 
decay or annihilation of dark matter particles. 

If hadronic production of mesons in the local ISM at energies well above TeV becomes 
the main source of high energy e± CRs near Earth, then the positron fraction should reach 
there the asymptotic value ~ 0.57 (Dado & Dar 2010; Cowsik et ah 2010) expected from the 
observed charge ratio /x + //x _ ~ 1.35 of high energy atmospheric and accelerator muons. 

The existence of a break near TeV in the combined e ± flux measured with H.E.S.S., 
still needs confirmation from independent and more precise measurements with instruments 
such as Fermi-LAT and the AMS. If verified, standard astroparticle physics may still explain 
it as shortly discussed in Section 6. Moreover, any alternative interpretation of its origin, 
including dark matter, must provide falsihable predictions to establish its validity. 
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Energy (GeV) 

Fig. 1.— Comparison between the high energy e + CR flux measured with AMS02 (Aguilar 
et al. 2014) and the secondary e + flux expected from hadronic interactions of the primary 
CR nucleons in the local ISM. 



















15 



Energy (GeV) 

Fig. 2.— Comparison between the flux of high energy e _ CRs measured with AMS02 (Aguilar 
et al. 2014) and the flux expected from Fermi acceleration in source of e _, s plus secondary 
production in/near source and in the ISM. 
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Fig. 3.— Comparison between the CR positron fraction near Earth as measured with AMS02 
(Accardo et al. 2014) and the positron fraction obtained from the predicted e ± CR fluxes 
shown in Figs. 1 and 2. 
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Fig. 4.— Comparison between the combined e 1 * 1 CR flux measured with AMS02 (Aguilar 
al. 2014) and the flux obtained from our standard particle astrophysics model. 
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Fig. 5.— The best fit ’minimal model’ positron fraction (Accardo et al. 2014) to the positron 
fraction measured with AMS02. 
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Fig. 6.— Comparison between the AMS minimal model (Accardo et al. 2014) best fit to the 
e + CR flux, Eq.(13), and the e + CR flux measured with AMS02 (Aguilar et ah 2014). 
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Fig. 7.— Comparison between the ’AMS minimal model’ (Accardo et al. 2014) best fit to 
the e _ CR flux, Eq.(14), and the e _ CR flux measured with AMS02 (Aguilar et ah 2014). 
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Fig. 8.— Comparison between the combined e 1 * 1 CR flux measured with AMS02 near Earth 
(Aguilar et al. 2014) and the flux obtained from the best fits of the ’AMS minimal model’ 
(Accardo et ah 2014) to their measured e + and e _ fluxes. 
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Fig. 9.— A power-law AE 317 fit to the combined high energy e 1 * 1 flux between 30-1000 
GeV measured near Earth with AMS02 (Aguilar et ah 2014). 
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Fig. 10.— A cutoff power-law fit to the combined high energy e 1 * 1 flux measured near 
Earth with AMS (full circles, Aguilar et al. 2014) and with H.E.S.S (squares, Aharonian et 
al. 2008,2009). The normalization of the H.E.S.S data was adjusted within their estimated 
systematic error to match the more precise AMS02 data below TeV (Aguilar et ah 2014). 




